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Abstract 

One  route  for  improving  the  performance  of  Li-based  batteries  is  to  optimize  the  carbon-based  electrode.  In  order  to  find  the  best 
carbon-based  materials,  the  specific  roles  of  the  molecular  and  solid-state  contributions  have  to  be  understood.  Here,  the  molecular 
contributions  are  analyzed.  A  semi-quantitative  method  is  proposed  to  compare  the  charge  storage  capability  of  polyaromatic  hydrocarbon 
molecules  (PAHs).  For  planar  PAHs,  the  ability  to  store  electrical  energy  is  found  to  be  to  a  large  extend  determined  by  a  single  parameter, 
that  is  the  electronic  hardness  (half  the  electronic  gap)  multiplied  the  number  of  carbon  atom  in  the  molecule.  A  compilation  of  results  for 
oligophenyls,  oligoacenes  and  medium-size  planar  systems  suggests  trends  in  the  dependence  of  the  energy  storage  capability  on  the  size 
and  shape  of  the  molecules. 
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1.  Introduction 

Lithium-ion  and  Li-metal  batteries  using  carbon-based 
electrodes  are  state-of-the-art  high  performance  power 
sources  [1-4],  The  low  standard  electrode  potential  of 
lithium  metal  («s— 3.05  V)  and  the  small  size  of  the 
lithium  ions  («s0.6A)  grant  high-energy  densities  and 
re-chargability  of  the  batteries;  thus  making  them  appropri¬ 
ate  for  portable  applications. 

One  route  to  improve  the  performance  of  Li-based  bat¬ 
teries  is  to  optimize  the  carbon  electrode  material  that  is  in¬ 
tercalated  with  lithium  ions  upon  charging  (Li-ion  battery) 
or  discharge  (Li-metal  battery)  of  the  device.  This  might 
be  achieved  by  a  better  understanding  of  the  fundamental 
mechanisms  and  the  energetics  underlying  the  intercalation 
process.  For  Li-metal  batteries,  in  particular,  the  two  key  pa¬ 
rameters  are:  firstly  the  amount  of  charge  stored  per  carbon 
atom;  secondly,  the  energy  gain  per  charge.  The  product  of 
both  determines  the  energy  gain  per  carbon  atom,  thus  the 
total  gravimetric  energy  density. 

High  levels  of  Li-intercalation,  up  to  about  LiC3.3, 
is  achieved  in  both,  amorphous  and  polycrystalline 
polyaromatic  hydrocarbon  (PAH)  materials  [5-9],  Q,o 
[10],  carbon  nanotubes  [11]  and  graphite  [12].  For 


*  Corresponding  author.  Tel.:  +46-13-281-217;  fax:  +46-13-288-969. 
E-mail  address:  raifr@ifm.liu.se  (R.  Friedlein). 

0378-7753/$  -  see  front  matter  ©  2003  Elsevier  B.V.  All  rights  reserved, 
doi:  10. 1016/j  .powsour. 2003.11.018 


poly(p-phenylene)-based,  as  well  as  polyacenic  semicon¬ 
ductor  (PAS)  materials  obtained  from  high-temperature 
treatment  of  polymers,  intercalation  levels  as  high  as  LiC~2 
are  reached  [6-13].  In  contrast  to  the  amount  of  Li-atoms 
intercalated,  the  amount  of  charge  transferred  is  usually 
hard  to  address.  Recent  quantum-chemical  calculations  sug¬ 
gest  an  incomplete  charge  transfer  for  the  LiC6  and  LiCg 
graphite  intercalation  compounds  [14,15]  as  well  as  for 
lithium  complexes  of  smaller  polyaromatic  hydrocarbons 
[16-19]  or  even  on  molecular  dimers  [20,21]. 

The  energy  stored  per  transferred  electron  is  best  de¬ 
scribed  by  the  formation  enthalpy  A H°  of  the  interca¬ 
lation  compound.  The  formation  enthalpy  of  an  isolated 
Li-molecule  complex  obtained  from  quantum  chemical 
calculations  does  not  match  directly  that  of  the  interca¬ 
lated  molecular  solid  [9,22,23].  An  accurate  estimate  of 
the  formation  enthalpy  of  the  Li-molecule  intercalation 
material  must  take  into  account  electronic  polarization  ef¬ 
fects  screening  the  electrostatic  interactions.  The  electronic 
polarization  characterizes  the  deformation  of  the  electron 
density  as  a  response  to  an  (external)  electric  field  [24].  One 
of  the  major  difficulties  lies  in  determining  the  Madelung 
energy  as  an  infinite  sum  of  polarization-screened  electro¬ 
static  interactions  between  positive  and  negative  charges. 
Recently  it  was  found  that  photoelectron  spectroscopy  gives 
many  information  on  the  lithium  intercalation  [9,23]:  the 
electronic  configuration  of  the  Li-atoms  in  organic  solids, 
the  formation  enthalpy  via  the  Born-Haber  cycle,  the  role 
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of  the  electronic  polarization  present  for  various  molec¬ 
ular  crystalline  phases,  and  even  a  first  estimate  of  the 
(polarization-screened)  Madelung  energy. 

In  the  present  paper,  we  identify  the  relevant  molecular  pa¬ 
rameters  for  the  electron-accepting  capability  of  molecules 
in  battery  applications.  Those  molecular  parameters  (elec¬ 
tron  affinity,  chemical  hardness)  are  normalized  to  the  num¬ 
ber  of  carbon  atoms  in  the  molecule,  or  alternatively  to  the 
mass  since  energy  densities  are  required  to  compare  various 
organic  solids. 


2.  Results  and  discussion 


2.1.  The  molecular  electron  affinity 


The  ability  of  a  molecule  to  accept  an  electron  is  deter¬ 
mined  by  the  electron  affinity.  The  adiabatic  electron  affin¬ 
ity  (EA)  of  a  free  molecule  is  defined  in  terms  of  the  to¬ 
tal  energy  difference  between  the  lowest  vibrational  state  of 
the  neutral  molecule  and  the  lowest  vibrational  state  of  the 
singly  negatively  charged  molecule.  The  vertical  electron 
affinity  is  the  energy  difference  between  the  neutral  and  neg¬ 
atively  charged  molecule  that  has  the  structure  of  the  neu¬ 
tral  molecule  in  the  ground  state.  Differences  between  the 
adiabatic  and  vertical  affinities  are  typically  of  the  order  of 
0. 1-0.2  eV  [25,26]. 

Density  functional  theory  (DFT)  provides  an  expression 
of  the  total  energy  of  a  molecule  as  a  function  of  the  number 
of  added  electrons.  The  total  energy  of  a  negatively  charged 
molecule  Etot  (x),  carrying  a  charge  x,  can  be  expressed  from 
the  total  energy  £tot(x  =  0)  of  the  neutral  molecule  via  a 
Taylor  series  [27-29]: 


Etot(N  =  No  -  x)  =  Eto\x  =  0)-x 


dEloi 


N=N0 


d2Etot 

8N 2 


R(x3) 


(l) 


N=N0 


where  R(x 3)  represents  the  higher-order  terms  of  the  series. 
Nq  stands  for  the  number  of  electrons  in  the  neutral  and  N  is 
the  charged  molecule.  The  higher-order  terms  are  expected 
to  be  small,  as  found  for  Cu-clusters  [29]  and  for  the  Cgo 
molecule  [30].  The  total  energy  of  the  charged  molecule  is 
then  a  quadratic  function  of  the  charge  state  x  as  shown  for 
Cgo-  Note  that  small  discontinuities  occur  only  after  filling 
of  shells.  Eq.  (1)  is  quantitatively  valid  for  x  between  —1  and 
1.  Hence,  its  use  for  higher  doping  levels  provides  only  a 
qualitative  tool  to  compare  PAHs  molecules.  For  non-integer 
charge,  Eq.  (1)  does  not  have  a  physical  meaning  for  isolated 
molecules,  although  it  can  be  used  for  open-systems  [31]  and 
eventually  molecular  solids  with  a  possible  partial  charge 
transfer  from  alkali  metal  atoms  to  the  molecular  Tr-system. 
The  advantage  of  using  this  formula  for  non-integer  charges 
is  to  provide  a  method  to  compare  the  charge  storage  ability 
of  various  molecules  for  a  certain  gravimetric  charge  density. 


In  DFT,  the  first  derivative  of  the  expansion, 
(dE/dN)N=N0,  considering  charge  addition,  is  defined  as 
the  chemical  potential  p.  of  the  electronic  system.  Like  in 
thermodynamics,  the  chemical  potential  expresses  how  the 
energy  changes  with  the  number  of  particles  in  the  system 
[32].  The  chemical  potential  can  be  seen  as  the  global  elec¬ 
tronegativity  for  a  molecule.  Indeed,  in  the  finite  difference 
approximation,  it  is  estimated  to  be  the  opposite  of  the 
average  between  the  ionization  potential  and  the  electron 
affinity  [29,33]: 

At*  — i(IP  +  EA),  r,  %  j  (IP  +  EA)  (2) 

which  is  nothing  but  the  Mulliken  electronegativity  for  atoms 
[34].  Half  of  the  second  derivative,  (1/2)  (32  Etot  /  dN2)  n=No> 
defines  the  absolute  hardness  i]  of  the  electronic  cloud  of 
a  molecule  (r/  is  always  positive)  [33],  The  nonchemical 
meaning  of  the  word  “hardness”  is  resistance  to  defor¬ 
mation  or  change.  In  DFT,  hardness  can  be  regarded  as 
the  resistance  of  the  chemical  potential  to  a  change  in 
the  number  of  electrons  [33].  It  is  also  related  to  the 
static  dipole  polarizability,  a  [35-37].  In  the  finite  differ¬ 
ence  approximation,  the  hardness  is  half  the  difference 
between  the  ionization  potential  and  the  electron  affinity 
(Eq.  (2)). 

Eq.  (2)  is  used  to  calculate  fi  and  i]  from  the  experimental 
values  of  IP  and  EA  for  a  number  of  polyaromatic  hydro¬ 
carbon  molecules  and  oligophenyls  (Table  1).  According  to 
Eq.  (1),  the  electron  affinity  of  a  molecule,  a  central  parame¬ 
ter  in  the  determination  of  formation  enthalpy  of  intercalated 
compounds  (Eq.  (2)),  can  be  written  in  terms  of  /x  and  rj\ 

EA  —  /j.  +  i]  (3) 

Since  the  chemical  potential  is  similar  (about  4  e  V)  for  most 
polyaromatic  hydrocarbon  molecules  [38,39]  (except  for  the 
fullerene  Cgo),  the  electron  affinity  variation  among  those 
molecules  is  directly  related  to  a  change  in  chemical  hard¬ 
ness.  A  significant  decrease  in  ri  is  observed  when  increasing 
the  size  of  the  molecular  Tr-system  (Table  1). 

According  to  Eq.  (2),  i]  appears  to  be  half  of  the  electronic 
gap  of  a  molecule.  Consequently,  the  larger  the  conjugated 
molecule,  the  smaller  the  electronic  gap,  the  smaller  the 
hardness  and  the  more  polarizable  the  electronic  cloud  [61]. 
As  expected,  larger  molecules  are  softer;  that  is,  it  is  easier 
to  charge  those  molecules  with  the  first  electron. 

2.2.  Optimal  electron-accepting  molecular  n-sy stems 

For  energy  storage  applications,  the  gravimetric  energy 
density  is  the  relevant  property  of  a  solid.  Among  the  poly¬ 
aromatic  hydrocarbons,  the  energy  storage  capability  of  each 
molecule  has  to  be  compared  for  the  same  gravimetric  charge 
density,  i.e.  for  the  same  number  of  charge  per  carbon  atom 
(neglecting  the  weight  of  the  hydrogen  atoms).  The  total 
energy  change  upon  charging  the  molecule  can  be  divided 
by  the  number  of  carbon  atoms,  Nq,  in  the  molecule.  Since 


Table  1 

Experimental  values  of  the  gas  phase  (vertical)  electron  affinity  and  (vertical)  ionization  potential  IPgas  of  some  polyaromatic  hydrocarbon  molecules  and  oligophenyls  obtained  from  electron  capture 
detector  (ECD)  or  photoelectron  spectroscopy  (PES)  measurements,  the  molecular  electronic  hardness  ij,  and  the  energy  cost  for  —0.1  charge  per  carbon  atom,  A£^ot  ( x/Nq  =  —0.1),  according  to  Eq.  (4) 


Molecule 

Electron  affinity  and  ionization  potential 

Chemical  potential  and  electronic  hardness 

Energy  cost  for 
— O.llel/C-atom 

Name 

Chemical  structure 

Nc 

EAgas  (eV)  IP«as  (eV) 

-n  (eV)  n  (eV)  Ncri  (eV) 

A£tot/A,c 

(x/Nc  =  -0.1)  (eV) 

Benzene 

Naphtalene 


Anthracene 


Tetracene 

(naphtacene) 

Pentacene 


Hexacene 

Pyrene 

Perylene 

Coronene 

Fullerene  C60 


Biphenyl 

p-terphenyl 


6 

-1.12  [40] 

9.24  [42] 

-1.15  [41] 

9.25  [43] 

10 

0.15  [44,45] 

8.15  [42] 

0.19  [40] 

0.20  [41] 

8.12  [43] 

14 

0.55  [25,44^16] 

7.40  [42] 

7.42  [47] 

7.47  [48] 

18 

1.15  [44,45] 

7.01  [42] 

1.03  [46] 

7.04  [48] 

22 

1.31  [46] 

6.64  [42] 

1.19  [49] 

6.61  [50] 

6.74  [48] 

26 

1.58  (theoretical 

6.44  (Ref.  [26];  some  results 

value  from  Ref.  [39]) 

from  other  work  cited  therin) 

16 

0.59  [44] 

7.41  (Ref.  [26];  some  results 
from  other  work  cited  therin) 

20 

0.98  [51,52] 

7.0  [53] 

26 

0.54  [49] 

7.25  (Ref.  [26];  some  results 

from  other  work  cited  therin) 

60 

2.8  [54] 

7.6  [55,56] 

2.65  [57] 

12 

0.13  [58] 

8.34  [59] 

18 

0.27  (theoretical 

7.9  (Ref.  [26];  some  results 

value  from  Ref.  [39]) 

from  other  work  cited  therin 

8.20  [60] 


4.06  ±  0.02 

5.19  ±  0.01 

31.1 

-0.09 

4.15  ±  0.10 

3.98  ±  0.02 

39.8 

-0.02 

3.99  ±  0.04 

3.44  ±  0.04 

48.2 

0.08 

4.06  ±  0.04 

2.97  ±  0.04 

53.5 

0.13 

3.97  ±  0.06 

2.72  ±  0.06 

59.8 

0.20 

*4.0 

62.6 

0.22 

«4.0 

«3.4 

54.4 

0.14 

*4.0 

**3.0 

60.0 

0.20 

«a3.9 

**3.4 

88.4 

0.49 

5.16  ±  0.07 

2.44  ±  0.04 

146.4 

0.95 

*4.2 

*4.1 

49.2 

0.07 

4.16  ±  0.10 

3.89  ±  0.08 

70.0 

0.28 
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higher  than  quadratic  terms  in  Eq.  (1)  are  negligible,  the 
energy  change  per  carbon  atom  is  then: 

A£tot(i)  _  Etot(N  =  N0  -  x)  -  Etot(N0) 

NC  Ac 

with  x/Nq  as  the  variable.  This  expression  provides  the  en¬ 
ergy  cost  per  carbon  atom  for  a  fixed  number  of  charges  per 
carbon  atom,  x/Nq-  Since  the  chemical  potential  is  similar 
for  planar  PAHs,  the  product  Net]  determines  the  ability  of 
the  molecule  to  store  energy. 

Table  1  displays  the  molecular  parameter  Ncq  and  the 
energy  cost  to  charge  the  molecules  with  0.1  electron  per 
carbon  atom.  At  this  charging  level,  energy  is  gained  by 
charging  the  benzene  and  naphtalene  molecules;  while  it 
costs  energy  for  the  other  molecules.  Molecular  systems 
with  small  Nc>]  require  less  energy  to  reach  the  same  degree 
of  charging.  Nc>]  appears  thus  as  the  relevant  parameter  to 
describe  the  energy  storage  capability  of  planar  poly  aromatic 
hydrocarbons;  and  it  is  plotted  in  Fig.  1  as  a  function  of  the 
size  of  the  molecule,  Nc- 

Several  trends  can  be  recognized: 

(1)  The  cost  to  charge  small  n-conjugated  systems  is  lower 
than  for  large  molecules  at  the  same  doping  level 
(xl  e  l/A'c).  Oligomers  are  favored  versus  polymers. 

(2)  Linearly  fused  aromatic  systems  (acenes)  need  less  en¬ 
ergy  to  become  charged  as  compared  to  non-linearly 
fused,  more  2D-like  molecules. 

(3)  Among  the  molecules  considered,  the  oligophenyls  have 
the  worst  electron-accepting  properties.  In  the  series 
from  benzene  to  p-terphenyl  Ac  >7  increases  linearly  with 
the  length  of  the  chain. 

(4)  In  contrast  to  the  excellent  ability  to  attract  a  single 
charge,  the  fullerene  C60  might  not  be  preferential  for 
energy  storage  applications,  where  the  amount  of  charge 
per  atom  matters. 


number  of  C-atoms  in  molecule  Nc 

Fig.  1.  The  product  of  the  molecular  electronic  hardness  r\  and  the  number 
of  carbon  atoms  in  the  molecule  Nc  of  acenes  (full  circles),  non-linear 
polyaromatic  hydrocarbon  molecules  (open  triangle)  and  oligophenyls 
(full  squares)  as  a  function  of  Nc  are  shown.  Lines  are  guides  for  the  eye. 


3.  Outlook 

A  semi-quantitative  method  is  proposed  to  compare  the 
energy  storage  ability  of  various  poly  aromatic  hydrocar¬ 
bon  molecules  at  low-charging  levels.  The  energy  storage 
in  molecules  is  discussed  in  terms  of  two  global  electronic 
properties  debited  in  density  functional  theory:  the  chemical 
potential  and  the  absolute  hardness.  Since  planar  PAHs  have 
similar  chemical  potentials,  the  absolute  hardness  multiplied 
by  the  number  of  carbon  atoms  in  the  molecule  appears  to 
be  the  relevant  molecular  parameter  to  determine  the  energy 
storage  in  molecules. 

Evidence  is  presented  that  small-  and  medium-size 
polyaromatic  and  phenyl-based  systems,  or  eventually  net¬ 
works  thereof,  might  have  the  highest  energy  storage  ability 
among  all  pure-carbon  systems.  This  statement  is  motivated 
by  the  lower  energy  cost  to  charge  those  molecules  for  a 
bxed  number  of  charge  per  carbon  atoms.  More  specif¬ 
ically,  for  a  given  (gravimetric)  charge  density,  linearly 
fused  aromatic  systems  (oligoacenes)  are  preferred  versus 
non-linearly  fused,  more  graphene-like  molecules  and  the 
oligo-  or  polyphenyls. 

To  be  able  to  extrapolate  safely  the  conclusions  obtained 
for  molecules  to  the  actual  Li-PAHs  intercalated  solids, 
Madelung  energy,  intermolecular  relaxation  and  electronic 
polarization  energies  should  also  be  taken  into  account;  this 
will  be  the  aim  of  a  future  work.  However,  the  observation 
that  small  Tr-conjugated  molecules  need  less  energy  to  be 
charged  at  a  debited  number  of  charge  per  carbon  atom  is  in 
agreement  with  their  high  level  of  energy  storage  measured 
experimentally  [9,24]. 
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